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Abstract. - The pseudogap Anderson impurity model provides a classic example of an 
essentially local quantum phase transition. Here we study its single-particle dynamics in the 
vicinity of the symmetric quantum critical point (QCP) separating generalized Fermi liquid 
and local moment phases, via the local moment approach. Both phases are shown to be 
characterized by a low-energy scale that vanishes at the QCP; and the universal scaling spectra, 
on all energy scales, are obtained analytically. The spectrum precisely at the QCP is also 
obtained; its form showing clearly the non- Fermi liquid, interacting nature of the fixed point. 



The celebrated Anderson impurity model (AIM) pQ, reviewed in [2], is the paradigm of 
quantum impurity physics. It describes a non-degenerate impurity with energy t\ and local 
interaction U, hybridizing via V\k = V to a non-interacting host with density of states p(u>) = 

5(u — £fc) (and u> — the Fermi level). The Hamiltonian in standard notation is 

H = ^2e k n ka + y^(e; + \Uh^ a )n- la + V ik (c\ a c k o + h.c). (1) 

h,(T o h.(T 

For a conventional metallic host, p(0) ^ 0, the problem is rather well understood the 
impurity spin being quenched and the system ubiquitously a Fermi liquid for all U. Its large- 
U behaviour is that of the Kondo effect, characterised by the single low-energy scale uk; as 
embodied famously in the Kondo resonance in the impurity single-particle spectrum D(u>), 
which scales universally in terms of uj/uk alone. 

The underlying physics is much richer if the host contains a power- law pseudogap 
p(u>) oc \o-i\ r , r > 0. This model has been studied extensively in recent years [3-12], especially 
for large-?/ where the low-energy subspace maps onto that of the Kondo model. It contains 
in general a non-trivial quantum phase transition at a finite U c = U c (r) ; the quantum critical 
point (QCP) separating a degenerate local moment (LM) phase arising for U > U c , from a 
'strong coupling' or generalized Fermi liquid (GFL) state [3-12] in which the impurity spin is 
locally quenched and a Kondo effect is manifest. As U — > U c — the Kondo scale u>k vanishes 
and the local spectrum again exhibits universal w/wK-scaling The QCP itself, where 

the weight of the Kondo resonance has just vanished, has been accessed very recently ^2] y i a 
a study of the local magnetization and susceptibility. The fixed point was found to be 
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interacting, exhibiting critical local moment fluctuations and an associated destruction of the 
Kondo effect very similar to those present at the local QCP of the Kondo lattice recently 
invoked to explain the behaviour of heavy fermion metals such as CeCug-ajAua; |13j . This 
provides topical motivation for further study of the model in the broad context of quantum 
phase transitions [14] . as too does its recently shown applicability to the issue of impurity 
moments in d- wave superconductors 

Here we study single-particle dynamics of the pseudogap AIM [6,8-11] via the local moment 
approach (LMA) |15|. a physically motivated many-body theory that has already led |8] 
to new predictions for the problem, borne out by direct numerical renormalization group 
(NRG) calculations [Jj]. Our aims here, where we focus on the large-?/ Kondo regime of 
the symmetric model, are threefold. First to obtain the universal scaling spectrum in the 
GFL phase: analytically and on all w/wK-scales, including but going far beyond the low- 
energy confines of traditional Fermi liquid theory. This is a non-trivial issue even for the 
metallic AIM (r = 0); where only recently has the scaling spectrum for uj/uik > 1 been 
shown |16j to be dominated not by power-law decays as long thought; but instead by long, 
slowly varying logarithmic tails. As we show, the situation is much richer in the pseudogap 
case, and significantly different from previous expectations [H1IH1- Second, in the LM phase 
where the Kondo scale = 0, we consider two questions: do the dynamics exhibit universal 
scaling close to the QCP on the LM side; if so what is the relevant low-energy scale, its 
physical origin and the form of the scaling spectra? Finally, we determine analytically the 
behaviour of the single-particle spectrum precisely at the QCP. Its low-energy behaviour is in 
turn shown to be intimately related to the high-energy asymptotics of the scaling spectra for 
both the GFL and LM phases; and its form shows clearly the non-Fermi liquid, interacting 
nature of the fixed point. 

Local moment approach. - We consider the symmetric QCP that arises in both particle- 
hole (p-h) symmetric and asymmetric models [SJEJ. To that end we here focus explicitly on 
the p-h symmetric AIM with e 4 = — -j (and local charge n; = ^ a {n- lc ,) = 1 V U); where 
the QCP separates GFL/LM phases for all < r < \ [5-9]. The host band is described 
by the (simplified) form p(u) — p \uj\ r 9(l — with bandwidth D = 1 (9(x): unit step 
function). We consider the T = impurity Green function G{uj) and hence local spectrum 
D{u) = — isgn(w)ImG(w); written conventionally as G(co) = [uj + — A(uj)~ £(w)] _1 , with lu + = 
u) + i0 + sgn(w) and £(o;) = S fl (w) — isgn(u)!] 1 (ui) the single self-energy (defined to exclude 
the trivial Hartree term). A(u) = Ar(w) — isgn(w)Ai(w) is the host-impurity hybridization, 
with Ai(w) = nV 2 p(uj) = A |w| r 6'(l — \lo\) where A = TrV 2 po; and Ar(w) follows by Hilbert 
transformation, A R (w) = -sgn(w)[/3(r)Ai(w) + 0(\w\)] with [3{r) = tan[|r]. 

The most illuminating expose of single-particle dynamics in the GFL phase resides in 
the modified spectral function !F{u}) = 7rAosec 2 ( : |r)|a;| r 'Z)(w) [8-10]; such that at the Fermi 
level J-(ui = 0) = 1 10 . The latter result is exact, generalizing to r > behaviour well 
known for the metallic AIM r = (where it amounts essentially trivially to satisfaction 
of the Friedel sum rule 0). It embodies the fact [E1E3 that the leading low-w behaviour, 
ttAqD(lo) ~ cos 2 ( : |r)|aj|~ r , is precisely that of the non-interacting limit: E R / J (w) vanishes as 
oj — > more rapidly than the hybridization (oc |w| r ), reflecting the perturbative continuity 
to the U = limit that in essence defines the GFL state. The Kondo resonance is directly 
apparent in !F(cu) [8-10] (see also Fig.l below). Determined by the Kondo scale u>k, it narrows 
with increasing U as the GFL/LM transition is approached and ojk 0. And the LMA 
predicts [Hj -^(w) itself to be universal in uj/ujk, as supported by NRG calculations 0- 

The LMA for the pseudogap AIM is detailed in (Hj. Rather than calculating S(w) directly 
it employs a two-self-energy description: the rotationally invariant G{uj) is written formally 
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as G(u) = \ Y,a G <y{^) with G a {u) = [uj+ - A(w) - E^u)]" 1 . The self-energies E a (oj) 
(= — E_ ff (— w) by p-h symmetry) are separated as Y,„(w) = — §J7|/i|+E (T (co'), into a static Fock 
piece with local moment that would alone survive at pure mean- field (MF) level; together 
with the dynamical E ff (u) = E^l^o-}] that naturally dominate the low-energy physics and 
are functionals of the broken symmetry MF propagators Q (T {u))= [ui + — A(oj) + cr-^- — 1 . In 
practice 8 , 15 : the LMA includes in E ct (oj) a non-perturbative class of diagrams that embody 
coupling of single-particle excitations to low-energy spin-flip processes, which is the essence of 
Kondo physics and provides a mechanism for dynamical communication between the broken 
symmetry MF states. These diagrams, shown e.g. in Fig 3 of |Hj, translate for E-[-(w) to 

/OO T 
— Imn + -(wi) e(ui)Ql (wi + u) + 8{-wi)Ql ( Wl + ^ ( 2 ) 
-oo 7T L 

where Q„{oS) denote the one-sided Hilbert transforms of the MF spectra 

D®(cu) = — —sgD.(u>)ImQ a (u>). II" 1 (oj) is the transverse spin polarization propagator, given at 

the simplest level employed here by an RPA-like p-h ladder sum (see e.g. Fig 3 of 

In describing the GFL phase the final, central idea behind the LMA is that of symmetry 
restoration |8l ll5lll6j : self-consistent restoration of the broken symmetry endemic at pure MF 
level, as embodied in Ej(0) = Ej(0) and hence (via p-h symmetry) 

E T (0) = E T (0)-i[/|/z|=0 (3) 

which ensures the Fermi liquid behaviour J-(0) = 1. This is achieved in practice joj for given 
U = U / A by varying the local moment from its pure MF value; which in turn introduces 
a low-energy scale to m = uj m (r) into the problem through a strong resonance in Imn -1 (u>) at 
u) = uj m . This is the Kondo or spin-flip scale, uj m = lo-^; and it sets the timescale r ~ h/ui m 
for symmetry restoration. In the LM phase by contrast, symmetry restoration Eq.(3) is not 
satisfied and w m = 0. ImlT^ (uS) here contains a pole at to = reflecting the zero energy 
spin-flip cost symptomatic of the locally degenerate phase; and the 'renormalized levels' E CT (0) 
are non-zero and sign-definite The GFL/LM phase boundary U — U c (r) may be obtained 
cither from the limit of solutions to Eq.(3) (as U — > U c — ); or, coming from the LM phase and 
yielding the same U c , from the condition that E CT (0) — > 0. GFL and LM phases are correctly 
found [5] to arise for all < r < |, and solely LM states for r > | and all U > 0. 

GFL scaling spectrum. Close to the QCP, !F(tS) scales in terms of Co = uj/u> m [HI 
E]. To obtain the scaling behaviour one formally considers finite uj in the limit u) m —* 0; 
thus projecting out non-universal features such as the Hubbard satellites. The hybridization 
A(w) = A(ojuj m ) then reduces to (Aquj^)' 1 A(uj) — — sgn(w)[/3(r) + i]|<D| r , and the 'bare' 
uj = ui m Q in G a (uj) = [u> — A(u) — Eg^cj)] -1 may be neglected (only r < | is relevant). Since 
J-(u) — 7rA sec 2 ( : |r)|a)| r 'wJ n -D(a;) exhibits scaling so too does uj^D(uj), given for ui > by 

A <^(c) = -Xi m ^[ (/3(r) + i)pr _ (AoO^S.H]- 1 . (4) 

Eq.(4) also shows that (Aow 1 ' n )~ 1 E cr (aj) must scale universally in terms of uj; to which we now 
turn within the LMA, following the approach of ^B] for the metallic AIM. Specifically we 
consider the strong coupling Kondo regime where 1 <C U (< U c (r) for the GFL phase). Since 
U c (f) ~ 7 as r — > 0, our analysis holds strictly for small r, although the central conclusions 
remain qualitatively valid save for r — > | (which the approach below does not handle). 

In strong coupling where the local moment saturates (|/i| — > 1), the transverse spin po- 
larization propagator entering Eq.(2) reduces to ilmIT + ~(a;) = S(ui — w m ). Hence E-|-(cj) = 
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U 2 Q7{u + uj m ), and in particular (S?j.(u;) =) E}(w) = 6(—[w + Lu m ])nU 2 D^lj + w m ). The 
leading low-w dependence of the MF D°(a;), which alone is relevant to the scaling behaviour, 
is D\(uo) ~ (4A /7rf/ 2 )|cj| r . Hence 

(A uj r J- 1 t}(oj)=9(-[cj + u; m })4\l + u,\ r (5) 

which indeed depends solely upon Co with no dependence on bare material parameters. The 
leading relevant low-o; behaviour of Re0r(w) is likewise given by |S] 



U 2 Reg7(uj) ~ A 



4 

— 7( r )i cj r 

irr 



(6) 



where 7(f) = 4/sin[7rr] ~ 4/irr ( and a uv-cutoff of order £) = 1 has been employed [S]). 

The dependence upon bare parameters of the Kondo scale w m (r) and critical U c (r) follows 
from symmetry restoration Eq.(3), viz (U 2 KeQ7 (w m ) =) S^(0) = in strong coupling. The 

critical {J c where u m = follows using Eq.(6) as U~ = which is exact as r — ► [HIEI - 
(For the corresponding Kondo model, the exchange coupling J relates to the AIM parameters 
under a Schricffcr- Wolff transformation by p^J = 8/wU ; whence the critical poJ c — 8/ttU c = r, 
as obtained originally by a large- N study of the pseudogap Kondo model [Hj). For U < U c , 
the [/-dependence of ui m follows from symmetry restoration as w m (r) = (I — [U /Uc}) 1 / 1 ' = 
(1 — ^C/) 1 / r . This recovers the exact exponential dependence uj m (r = 0) = exp(— jU) for 
the r = metallic AIM 2 , in turn showing the 1/r exponent of u> m (r) to be exact as r — > 0. 
Finally, we require the (Independence of S^w) = — S?(0), given using Eq.(6) by 

(AoO- 1 ^^) = _ 7 (r)[|l + cr - 1] (7) 

which as required scales solely in terms of Q. Eqs.(5,7) for E-f(w) (= — S^(— u>)), together with 
Eq.(4), determine the entire LMA scaling spectrum; given for Co > (CF(u>) = CF(—u)) by 



1 ~ 2r 

^(w)cos 2 [fr] = — — — - 2 + 



(8) 

4^-1)11 -in 



\P(r)& - 7(r)(|w - l| r - l)] 2 + w 2 »-[l + 40(£ - - i|f 

with /3(r) = tan^r] ~ ?r and 7 (r) ~ — . Eq.(8) reproduces the scaling spectrum obtained 
numerically in [8j, fully quantitatively as r — ► 0. For r = 0,0.05,0.1 and 0.15 it is illustrated 
in Fig.l on all relevant w-scales, now discussed. The leading Co — > behaviour is J~{oj) ~ 
1 — ArCo 1 ^ 1 ' — (4/ 7r) 2 a> 2 ^ 1— r - ) , with an |w| 1_r cusp (Fig.l inset) whose weight vanishes as r — > 0; 
recovering for r = the conventional Fermi liquid form for the metallic AIM, J-(oS) — 1 ~ Co 2 . 

We turn now to the spectral 'tails', arising formally for d 3> 1 ^7]. Based on numerical 
fitting of LMA and NRG results it was argued in 9| that for all < r < h , as hitherto believed 
for r = ( JH| an d references therein), the power-law decay uj^DIlo) ~ |a)| -1 / 2 arose. For the 
metallic AIM however, this is now known to be wrong: recent work |16| shows the tail decay 
to be far slower, in fact logarithmically so for all |cD| > 1. And the situation is more subtle for 
r / 0. For | tZ> | 3> 1, but for \Co\ r = c rln l"l otherwise arbitrary, Eq.(8) gives 

1 f 1 



HUJ) ~ 2 \[£(i_e-rHa|)]2 + i + [A(i_ e -HnM )]2 + 2 5 / 
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Fig. 1 - LMA scaling spectrum for the GFL phase (Eq.(8)): T(lo) vs. u> — Lu/um (oj m = i^k) on all 
relevant energy scales; for r — (dotted), 0.05 (dashed), 0.1 (long-dashed) and 0.15 (solid). 



(where the leading low-r behaviour is used). For r = 0, where -(1 — e _rln '"') = m|<D|, this 
reduces to the known result (Eq.(4.2) of ^5]) which in practice agrees quantitatively with NRG 
calculations for ii > 5 ^Hl- That same form clearly obtains for r / provided 1 -C \u>\ <C e 1 / 1 "; 
and for r < 1 'log-tails' thus appear to dominate the spectrum, see Fig.l. But for r ^ the 
ultimate asymptote of approached for \uj\ t 3> 1 (see Fig.l), is a non-zero constant: 

~ (37r 2 /16)r 2 = (3tt 2 /16)[po^c] 2 ; and hence 

7rA <£>M ~ ^r 2 \Cu\- r (10) 

- where, since u = ui/u> m , we note that the Kondo scale u> m 'drops out' of the tail. Cor- 
responding results for the single self-energy £(w) are also readily obtained via S(w) = u> — 
A(w) — 1/G(w); such that (A o;J n )~ 1 I](a;) scales universally in Q. For 1 <C \u>\ <C e 1 / 1 " one 
finds the intermediate behaviour (AoW^ 1 ) _1 S I (o;) ~ (16/37r 2 )ln 2 (|<D|), just as for r = 
but its ultimate asymptotic form for \uj\ r ^> 1 is 

(AqO- 1 ^^) ~ (16/37r 2 r 2 )|^r (11) 

(with S r (cj) ~ —sgn.(oj)/3(r)Ti I (u)), from which w m again drops out. 

LM scaling spectrum. - But what of the LM phase U > U c (r), throughout which (as 
above, see also |S]) the Kondo spin-flip scale uj m = 0; and where for \ui\ — > 0, D(ui) is 
known to vanish as D(u>) oc \lu\ t 6,8,9]? Here, since oj m — 0, S-f(w) = U 2 Q^(u>) in strong 
coupling (U ~> 1), and hence A^" 1 E|(w) = &{— uj)4\uj\ r {i.e. Eq.(5) with uj m = 0); likewise 
Aq^E^cj) - tf (0)] = -7(r)|w| r (Eq.(7) with w m = 0). But in the LM phase symmetry 
is not restored, hence the 'renormalized level' S^(0) = — ^ + U 2 RcQ~^ (0) ^ 0; and is given 
via Eq.(6) by Aq £?(0) = — |[f/ — U c ], vanishing linearly as U — > t/ c +. It is now simple to 
show that the LM phase is itself characterized by a non-zero low-energy scale ui\ = u>i(r) that 
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Fig. 2 - Schematic in the (|u;|, U) plane of spectral behaviour in the vicinity of the QCP, with ui = 
the Fermi level. The crossover scale oj* (see text) is shown as a dashed line. 



is determined by the renormalized level and given by lo\ = [(A 7(r)) 1 |S^(0)|] 1 / r ; whence 

u>\ = {[U /U c ] — l) 1 ^ using U c (r)/2j(r) = 1 for small r (which form, for A < r < 1, has also 
been found in a recent NRG study close to the asymmetric QCP of the model). Close to 
the (symmetric) QCP we find that uj[D(u)) indeed scales universally in terms of uj — u>/u>i, 
and is given explicitly by: 

TrAowTDfw) = ^-r 2 , - „ (12) 

1 V ; i 6 ( 1 + |£|r)2 V ) 

The LM and GFL phases are thus each characterized by a low-energy scale a;* = | [U/U c ] — 1 1 1 / r 
for small r, reducing to u m [uj\] for U < U c [> U c ] in the GFL [LM] phase; and vanishing as 
the QCP is approached from either phase. Eq.(12) recovers the known lowest-w behaviour 
nA D(u)) ~ 3(Ao/S| i (0)) 2 |tj| r . And for \uj\ r ^> 1, noting that u m = = wj as above, 
Eq.(12) for the LM phase reduces precisely to Eq.(lO) for the GFL phase: the scaling spectra 
of the two phases share common tails. 

Quantum critical behaviour. - Precisely at the QCP, the low-energy scale = 0. In 
consequence the low-w behaviour of D(ui) is given by 

nA D(cj; U = U c (r)) = ^r 2 M" r = ^(po JcfH^ (13) 

as follows by taking the limit = in either Eq. (10) or (12); or equivalently by employing 
= from the outset. The important point is that at the QCP Eq.(13) holds as \uj\ — > 0, 
with a prefactor to the \uj\~ r divergence that is renormalized from its non-interacting value 
(of cos 2 [ : |r] ~ 1) to precisely the high-energy plateau value of in the GFL phase, found 
above. Physically this embodies the fact that as U — > U c {r), the Kondo resonance in !F{lS) 
collapses at the QCP, leaving only the tails and producing a featureless F{uj) that is no longer 
pinned to unity at the Fermi level, ^(0) = (37r 2 /16)r 2 . In this way (and hence D(u})) at 

the QCP is distinct from that approaching it, where u> m = is small but strictly non-zero and 
any rescaling in terms of it recovers Eq.(8) (and hence .F(O) = 1). The non-[G]FL, interacting 
nature of the fixed point ^21 is clearly evident; embodied both in .F(O) ^ 1 and, relatedly, in 
the -> behaviour (Ao)" 1 !] 1 ^) = (16/37r 2 r 2 )|w| r (Eq.(ll) with oj m = 0), whose decay 
to zero with the same power as the hybridization vitiates the necessary condition 8 , lO] for 
adiabatic continuity to the non-interacting limit and hence a GFL state. 

In summary we have considered a local moment approach to the pseudogap Anderson 
model close to the symmetric quantum critical point; focusing on T = single-particle dy- 
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namics, in which the collapse of the Kondo resonance as the QCP is approached is directly 
manifest (the present approach does not enable us to handle on a comparable footing the local 
magnetic susceptibility x['oc considered in H2J). Local impurity spectra close to the QCP have 
been shown to exhibit scaling in terms of a low-energy scale, a;* = |[t//f7 c ] ~ l| 1/,r for small r, 
that vanishes as the QCP is approached from either phase; and scaling spectra in both phases 
have been obtained analytically. The situation is summarised in Fig. 2, representing a natu- 
ral complement to the conventional picture in the (T,£/)-plane that schematises the crossover 
behaviour of static properties in the vicinity of a QCP ^31- Here by contrast, crossovers in 
dynamics are exemplified in the (|ct>|,J7)-plane at T = 0: for \u>\ <C u>* the low-energy spectral 
characteristics of the GFL and LM phases are distinct, while for \u>\ 3> each reduces to the 
common form ttAqD^) ~ (37r/16)(po>/c) 2 M _r ; such that the latter behaviour, precisely that 
of the QCP itself, ever increasingly dominates the dynamics of either phase as the transition 
is approached and vanishes. The pseudogap AIM is itself a paradigm for understanding 
generic aspects of local quantum phase transitions, see e.g. [11-13]; and we are currently ex- 
tending our approach to encompass the asymmetric case and the influence of a local magnetic 
field. 

* * * 
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